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Abstract 

Background: Soybean is an important crop that provides valuable proteins and oils for human use. Because 
soybean growth and development is extremely sensitive to water deficit, quality and crop yields are severely 
impacted by drought stress. In the face of limited water resources, drought-responsive genes are therefore of 
interest. Identification and analysis of dehydration- and rehydration-inducible differentially expressed genes (DEGs) 
would not only aid elucidation of molecular mechanisms of stress response, but also enable improvement of crop 
stress tolerance via gene transfer. Using Digital Gene Expression Tag profiling (DGE), a new technique based on 
lllumina sequencing, we analyzed expression profiles between two soybean genotypes to identify drought- 
responsive genes. 

Results: Two soybean genotypes — drought-tolerant Jindou21 and drought-sensitive Zhongdou33 — were subjected 
to dehydration and rehydration conditions. For analysis of DEGs under dehydration conditions, 20 cDNA libraries 
were generated from roots and leaves at two different time points under well-watered and dehydration conditions. 
We also generated eight libraries for analysis under rehydration conditions. Sequencing of the 28 libraries produced 
25,000-33,000 unambiguous tags, which were mapped to reference sequences for annotation of expressed genes. 
Many genes exhibited significant expression differences among the libraries. DEGs in the drought-tolerant genotype 
were identified by comparison of DEGs among treatments and genotypes. In Jindou21, 518 and 614 genes were 
differentially expressed under dehydration in leaves and roots, respectively, with 24 identified both in leaves and 
roots. The main functional categories enriched in these DEGs were metabolic process, response to stresses, plant 
hormone signal transduction, protein processing, and plant-pathogen interaction pathway; the associated genes 
primarily encoded transcription factors, protein kinases, and other regulatory proteins. The seven most significantly 
expressed (|log 2 ratio| > 8) genes— Glyma15g03920, Glyma05g02470, Glyma15g15010, Glyma05g09070, 
Glyma06g35630, Glyma08g 12590, and Glymal 1 gl 6000 — are more likely to determine drought stress tolerance. The 
expression patterns of eight randomly-selected genes were confirmed by quantitative RT-PCR; the results of 
QRT-PCR analysis agreed with transcriptional profile data for 96 out of 128 (75%) data points. 
(Continued on next page) 



* Correspondence: zhouocri@sina.com; wenbinli@yahoo.com 

^il Crops Research Institute of Chinese Academy of Agriculture Sciences, 

Wuhan 430062, China 

3 Key Laboratory of Soybean Biology in the Chinese Ministry of Education, 

Northeast Agricultural University, Harbin 150030, China 

Full list of author information is available at the end of the article 

O© 2013 Chen et al.; licensee BioMed Central Ltd. This is an open access article distributed under the terms of the Creative 
BiolVlGCl C6ntTcll Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



Chen et al. BMC Genomics 2013, 14:687 
http://www.biomedcentral.eom/1 471 -21 64/1 4/687 



Page 2 of 19 



(Continued from previous page) 

Conclusions: Many soybean genes were differentially expressed between drought-tolerant and drought-sensitive 
genotypes. Based on GO functional annotation and pathway enrichment analysis, some of these genes encoded 
transcription factors, protein kinases, and other regulatory proteins. The seven most significant DEGs are candidates 
for improving soybean drought tolerance. These findings will be helpful for analysis and elucidation of molecular 
mechanisms of drought tolerance; they also provide a basis for cultivating new varieties of drought-tolerant 
soybean. 

Keywords: Soybean, Dehydration, Digital gene expression tag profile, Rehydration, Differentially expressed genes, 
Quantitative RT-PCR, Transcription factors, Protein kinases, Regulatory proteins 



Background 

Drought stress is one of the major environmental factors 
adversely influencing plant growth and development [1]. 
In the face of globally-increasing areas of aridity and 
semi-aridity, the improvement of plant drought toler- 
ance is crucial [2], The primary plant response to water 
deficit is reduction in water content, leading to stomatal 
closure and/or decreased transpiration. Because stomatal 
closure in turn leads to reduction in chlorophyll content 
and C0 2 assimilation/photosynthetic rates, stomatal sta- 
tus is currently a topic of intense research interest. Early 
symptoms of water deficit in plants include leaf wilting 
and rolling and leaf etiolation. Plants adapt to drought 
by means of various processes, including stress stimulus 
response, signal perception, signal transduction, and 
stress-responsive gene expression, with many morpho- 
logical and physiological changes at molecular and cellu- 
lar levels providing protection against biotic and abiotic 
stress damage [3]. 

Much biological and genetic diversity exists among 
different plant species with respect to adaptation to 
drought stress. Many of these stress adaptive mecha- 
nisms are not completely understood. Identification of 
novel drought-responsive genes from drought-tolerant 
plants and elucidation of their roles in drought adapta- 
tion will help improve the tolerance of drought-sensitive 
plants [4,5]. 

Drought tolerance is a complex quantitative trait, with 
many genes induced in response to drought. Drought- 
responsive genes were initially identified in Arabidopsis 
using microarray analysis. The products of these genes 
can be divided into two main classes [6]: 1) functional 
proteins, such as water channel proteins, detoxification 
enzymes, late embryogenesis abundant (LEA) proteins, 
key enzymes for osmolyte biosynthesis (praline, sugars), 
and proteases [7-9], and 2) regulatory proteins, including 
transcription factors (TFs), protein kinases (PKs), and 
proteins associated with phospholipid metabolism and 
abscisic acid (ABA) biosynthesis. 

Production of DNA sequencing reads using the 
Illumina/Solexa platform is currently widespread [10]. 
Illumina/Solexa sequencing-based Digital Gene Expression 



Tag profile (DGE) technology can be used to comprehen- 
sively, economically, and rapidly detect changes in gene ex- 
pression [11]. Because of its high quantitative accuracy, 
repeatability, and wide detection threshold [12], DGE tech- 
nology has been widely used to study differentially 
expressed genes (DEGs) in various fields of biology, includ- 
ing human disease, animal immunization, and plant re- 
search [13-15]. 

Soybean (Glycine max [L.] Merr.) is an important 
crop plant because of its seed protein and oil content. 
Soybean quality and yield are adversely influenced by 
drought stress [16,17], which can cause yield losses of 
approximately 40%. To overcome these limitations, 
the identification of key drought-responsive genes is 
needed. With the completion of soybean genome se- 
quencing, the soybean genome has been well character- 
ized [18]. Multiple genes have been annotated using 
functional genomic methods and their functions tenta- 
tively identified. 

In this study, a large number of drought-responsive 
genes differentially expressed between drought-tolerant 
and drought-sensitive genotypes were identified using an 
Illumina/Solexa sequencing system. Some of the DEGs 
uncovered in this study should aid efforts to understand 
soybean drought tolerance mechanisms. 

Methods 

Preliminary screening of soybean materials under 
dehydration 

To obtain drought-tolerant and drought-sensitive mate- 
rials, 57 soybean genotypes were collected from different 
areas of China (Additional file 1). Sterilized seeds were 
germinated on moistened filter paper for 5-6 d at 28°C 
and 60% humidity under a 16 h/8 h (light/dark) cycle. 
The seedlings were then transferred into plastic boxes 
and cultured hydroponically. Drought treatment was 
carried out as follows: when the first trifoliolate leaves 
unfolded, plants to be treated were transferred into plas- 
tic boxes without water. They were allowed to dry for 0, 
0.5, 1, 2, 4, 6, 8, 10, and 12 h at 28°C at 60% humidity 
under a 16 h/8 h photoperiod. Control plants were 
maintained in water under the same conditions. Drought 
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symptoms (leaf wilting and rolling in the top of the can- 
opy) could be readily observed. Wilting was rated on a 
scale of 0 to 100 (0 = none, 20-40 = slight, 40-60 = mod- 
erate, 60-80 = severe, and 100 = plant death) [19]. Further 
screening was conducted using the related drought- 
tolerance indices of relative water content (RWC) and 
relative electrical conductivity (REC). RWC and REC 
measurements were performed according to published 
methods [20,21]. 

Sample preparation and library construction 

Drought-tolerant Jindou21 and drought-sensitive Zhongdou33 
were selected for sequencing. Seeds were germinated 
on filter paper for 5-6 d. Seedlings were transferred to 
large plastic plates filled with water and grown under 
greenhouse conditions (28°C, 16 h/8 h photoperiod, 
80 (imol m" 2 s" 1 photon flux density, and 60% relative 
humidity). Root and leaf tissues were collected separ- 
ately when the first trifoliolate soybean leaves unfolded. 
Drought treatment was carried out as follows [22]: 
plants were transferred onto filter paper to soak up 
water, and then allowed to dry for 0, 2, and 10 h at 
28°C under a 16 h/8 h photoperiod, 13 [irnol m" 2 s" 1 
photon flux light intensity, and 60% relative humidity. 
Control plants were maintained in water for 0, 2, and 
10 h under the same conditions (Figure 1). In addition, 
seedlings that had been allowed to dehydrate for 2 h 
were rehydrated separately for 0.5 h and 2 h. Root 
and leaf tissues of dehydrated, rehydrated, and control 



plants were separately collected, with three biological 
replicates, for sequencing. 

Twenty-eight cDNA libraries were generated for se- 
quencing and expression profile analysis: 20 were 
constructed from dehydrated plants, and 8 were gener- 
ated from plants undergoing rehydration after a 2-h de- 
hydration treatment (Table 1). 

Illumina/Solexa sequencing and clean tag library 
formation 

Sequencing and library formation were performed using 
an Illumina Gene Expression Sample Prep Kit and a 
Solexa Sequencing Chip (flow cell), with primary instru- 
mentation consisting of an Illumina Cluster Station and 
an Illumina HiSeq 2000 System. Raw sequences, includ- 
ing 3' adaptor fragments, low-quality sequences, and 
several types of impurities, were generated as detailed in 
Additional file 2. Raw sequences were transformed into 
clean tags by (1) trimming 3' adaptor sequences from 
the 49-nt raw reads, and then removing (2) empty reads 
(reads containing 3' adaptor sequences without tags), (3) 
low quality tags (tags with unknown sequences 'N'), (4) 
tags with lengths of other than 21 nt, and (5) tags with a 
copy number of one (probably due to sequencing error). 
In the remainder of this paper, the term 'total clean tags' 
corresponds to the number of clean tags generated, 
while 'distinct clean tags' refers to the number of clean 
tag types produced. Tag preparation principles and steps 
are further detailed in Additional file 2. 
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Table 1 28 libraries for sequencing under dehydration and rehydration conditions 



Process 



Library 





JWOL 


JWOR 


JD2L 


JW2L 




Jindou21 -water 


Jindou21 -water 


Jindou21 -dehydration 


Jindou21 -water 




0 h-leaves 


0 h-roots 


2 h-leaves 


2 h-leaves 


P)ph\/fl ration 

L/CI 1 VUI CI LI Ul 1 


JD2R 


JW2R 


JD10L 


JW10L 




JIIIUUUZI utriiyu I d UUI I 


JIIIUUUZI Vvdlcl 


jiiiuuuzi utri iyu 1 d LIUI I 


li nHni i ") 1 ->A/2tor 
JIIIUUUZI VVdLcl 




Z II IDOLS 


Z II lUULS 


1 U 1 1 Itrdvtrb 


1 U 1 1 Itrdvtrb 




JD10R 


JW10R 


ZW0L 


ZWOR 




Jindou21- dehydration 


Jindou21 -water 


Zhongdou33-water 


Zhongdou 33- water 




1 0 h-roots 


1 0 h-roots 


0 h-leaves 


0 h-roots 




ZD2L 


ZW2L 


ZD2R 


ZW2R 




Zhongdou33- 


Zhongdou33-water 


Zhongdou33- 


Zhongdou33-water 




HohvH r ration 
Utrl lyUldUUI 1 










2 h-leaves 


2 h-leaves 


dehydration 2 h-roots 


2 h-roots 




ZD10L 


ZW10L 


ZD10R 


ZW10R 




Zhongdou33- 


Zhongdou33-water 10 h- 


Zhongdou33- 


Zhongdou33-water 10 h- 




dehydration 


leaves 


dehydration 


roots 




10 h-leaves 




1 0 h-roots 






JR0.5 L 


JR0.5R 


JR2L 


JR2R 




Jindou21- rehydration 


Jindou21- rehydration 


Jindou21- rehydration 


Jindou21- rehydration 




0.5 h-leaves 


0.5 h-roots 


2 h-leaves 


2 h-roots 


Rehydration after 2-h 


ZR0.5 L 


ZR0.5R 


ZR2L 


ZR2R 


dehydration 


Zhongdou33- 


Zhongdou33-rehydration 


Zhongdou33- 


Zhongdou33-rehydration 




rehydration 




rehyd ration 






0.5 h-leaves 


0.5 h-roots 


2 h-leaves 


2 h-roots 



Each cell lists genotype (Jindou21 or Zhongdou33), treatment (water control, dehydration, or rehydration), treatment duration, and tissue sampled. 



Gene expression annotation 

A virtual tag library containing all 17-nt + CATG se- 
quences was generated using the SoyBase soybean gen- 
omics database (http://www.soybase.org/). All clean 
tags were mapped to these reference sequences, with 
only tags having mismatches < 1 bp retained. After fil- 
tering out clean tags mapping to reference sequences 
from multiple genes, the remaining clean tags were 
designated as unambiguous clean tags. The number of 
unambiguous clean tags for each gene was calculated 
and normalized to the number of transcripts per mil- 
lion clean tags (TPM). 

RNA isolation, DNase I treatment, and cDNA synthesis 

Total RNA was isolated using Trizol reagent (Invitrogen, 
USA) and quantified with an Epoch Multi-Volume Spec- 
trophotometer system. Each 6-|ig total RNA sample was 
digested with RNase-free DNase I (Fermentas, Canada) 
to remove genomic DNA contamination. After DNase I 
treatment, RNA concentration was checked again; 
A260/A280 values of all RNA samples used in this study 
ranged from 1.8 to 2.2 and A260/A230 ratios of all sam- 
ples were above 2.0. First-strand cDNA synthesis was 



performed in 20- (il reaction volumes using a DNA Syn- 
thesis kit (Promega, USA). 

Housekeeping gene selection and quantitative real-time 
PCR analysis 

Housekeeping genes are used as reference genes to 
quantify gene expression. Although reference genes 
should ideally be expressed at a stable level throughout 
the plant and not be influenced by exogenous treatments 
[23,24], studies have found that transcription levels of 
some housekeeping genes vary considerably in response 
to changes in experimental conditions and across differ- 
ent tissue types [25,26]. For this reason, we selected 
five housekeeping genes to evaluate levels of gene ex- 
pression: TUA5, ACT2/7, UBQ10, UKN2, and HDC 
(Additional file 3). Stability of reference gene expression 
was analyzed using geNorm software (v3.50). We created 
a ranking of stability measure values by stepwise exclu- 
sion of the least stable housekeeping gene, and then 
used the relatively most stable housekeeping gene to 
normalize expression levels of selected genes [27]. 

Primers were designed using Primer 5.0 software. Ex- 
pression levels for all candidate genes were computed 
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based on the stable expression level of the reference 
gene. Quantitative real-time PCR (QRT-PCR) was 
performed in 96-well plates on a Bio-Rad iQ5 Real-Time 
PCR Detection system using SuperReal PreMix (SYBR 
Green) reagents (Tiangen, China). Reactions were car- 
ried out in 10- ul volumes containing 0.3 uM final con- 
centrations of each primer. The QRT-PCR thermal 
profile consisted of 95°C for 15 min, followed by 40 cy- 
cles of 95°C for 10 s and 60°C for 32 s. Dissociation 
curves were obtained from a thermal melting profile 
generated under a final PCR cycle of 95°C for 30 s 
followed by a constant increase in temperature from 60°C 
to 95°C. Threshold values were empirically determined 
based on the observed linear amplification phase of all 
primer sets. Sample cycle threshold (ct) values were 
standardized for each template based on reference gene 
control primer reaction, and the 2 _AACT method was 
used to analyze relative changes in gene expression [28]. 
Three replicate reactions per sample were used to ensure 
statistical credibility. 

Results 

Screening of soybean genotypes for drought tolerance 

To identify soybean genotypes with drought tolerance, 
57 genotypes were screened at the seedling stage under 
water-deficit conditions (Additional file 1). Based on 
preliminary screening experiments, we selected two 
drought- tolerant genotypes (Jindou21 and Tianlongl) 
and two drought-sensitive genotypes (Zhongdou33 and 
Zhongdou34). Further screening experiments revealed 
significant differences in RWCs and RECs among the 
four soybean genotypes under dehydration and rehydra- 
tion conditions. RWCs of all four genotypes declined 
during dehydration and increased after rehydration, but 
the RWC of Jindou21 was always higher than that of the 
other genotypes under drought stress and quickly recov- 
ered to control levels after 2 h of rehydration (Figure 2). 





Jindou 21 




Zhongdou 33 




Zhongdou 34 




Tianiong No. 1 




0 D2 R-D2 D10 R-D10 

Figure 2 Relative water content of four genotypes under 
dehydration for 0, 2, and 10 h and under rehydration. D2 = 2 h 

dehydration; R-D2 = rehydration after 2 h dehydration; D10= 10 h 
dehydration; R-D10 = rehydration after 10 h dehydration. 



REC measurements indicated that RECs of all four 
genotypes increased under dehydration. The extent of 
increase was lowest in Jindou21 and highest in 
Zhongdou33, which had the poorest recovery after 2 h 
dehydration. RECs of all four genotypes were elevated 
during rehydration following 10 h dehydration (Figure 3), 
indicating that prolonged dehydration may have dam- 
aged the seedling leaves. Based on these results, we 
choose the drought- tolerant genotype Jindou21 and the 
drought-sensitive genotype Zhongdou33 for DGE. 

Illumina/Solexa sequencing evaluation 

In this study, an Illumina/Solexa system was used for 
tag sequencing analysis. Twenty-eight samples from 
Jindou21 and Zhongdou33 were prepared for sequen- 
cing: 20 samples were obtained from 0, 2, and 10 h de- 
hydration treatments and 8 samples were derived from 
0.5 h and 2 h rehydration treatments. After identifying 
DEGs between control (water) and dehydration treat- 
ments, DEGs between the two genotypes were gener- 
ated. The number of raw sequence tags generated from 
each of the 28 libraries ranged from 5.8 to 6.2 million, 
and the number of raw tags producing distinct sequences 
ranged from 0.3 to 0.8 million (Additional file 4). To 
assess Illumina/Solexa sequence quality, a sequencing 
saturation analysis was performed to check whether the 
number of detected genes increased with increasing 
amounts of sequence data. After generation of about 
three million tags, the number of detected genes barely 
continued to increase (Additional file 5); at this point, 
the library was saturated and contained sufficient infor- 
mation for gene expression analysis. 

After trimming and filtering the raw sequence data, 
there were many distinct tags generated from the 28 li- 
braries. The distribution of the various tag abundance 
categories between total and distinct tag counts was 
consistent across all libraries. Matching tags to genes is 
necessary for sequence prediction and determination of 
molecular mechanisms underlying gene expression. 
When clean tags were BLAST-searched against available 
genome data (http://phytozome.net/), more than 75% 
mapped to known genes. The number of tag-mapped 
genes represented in each library ranged from 30,155 to 
29,706 (Additional file 4). 

Identification of DEGs between the two treatments 

To judge the significance of differences in expressed 
genes, we used three criteria: FDR < 0.001, |log 2 ratio | > 1, 
and P- value < 0.01. In this study, many genes were differ- 
entially expressed between control (water) and dehydra- 
tion conditions, and between dehydration and rehydration 
treatments. For example, in the drought-tolerant genotype 
Jindou21, there were 3,337 DEGs between (JW2LnJW0L) 
vs. JD2L, 4,968 DEGs between (JW2RnJW0R) vs. JD2R, 
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12,340 DEGs between (JWlOLnJWOL) vs. JD10L, and 
6,449 DEGs between (JWlORnJWOR) vs. JD10R. In the 
drought-sensitive genotype Zhongdou33, 4,021 DEGs 
were detected between (ZW2LnZW0L) vs. ZD2L, with 
3,460 DEGs observed between (ZW2RnZW0R) vs. ZD2R, 
12,137 DEGs between (ZWIOLnZWOL) vs. ZD10L, and 
8,857 DEGs between (ZWIORnZWOR) vs. ZD10R. There 
were also many DEGs between the two treatments. In the 
drought-tolerant genotype, there were 5,191 DEGs ob- 
served between JR0.5 L vs. JD2L, 6,578 DEGs between 
JR0.5R vs. JD2R, 5,552 DEGs between JR2L vs. JD2L, and 
6,100 DEGs between JR2R vs. JD2R; 2,286 DEGs were 
found between ZR0.5 L vs. ZD2L, 3,968 DEGs between 
ZR0.5R vs. ZD2R, 6,665 DEGs between ZR2L vs. ZD2L, 
and 5,483 DEGs between ZR2R vs. ZD2R in the drought- 
sensitive genotype (Figure 4). A Venn diagram illustrating 
the division of DEGs into different groups according to 
dehydration time, expression site, and rehydration process 
is shown in Figure 5. 

Further analysis of DEGs between the two genotypes 

Based on results observed between treatments, differ- 
ences between the two genotypes were further compared 
with a focus on DEGs found only in the drought-tolerant 
genotype. Under dehydration, a number of genes were 
expressed only in the drought-tolerant genotype: 518 in 
leaves (Additional file 6), 614 in roots (Additional file 7), 
and 24 in both leaves and roots. Although some of these 
genes could be assigned putative functions, the functions of 
most were unknown. In Jindou21 under rehydration condi- 
tions, 563 genes were expressed only in leaves (Additional 
file 8), 1,720 in roots (Additional file 9), and 65 both in 
leaves and roots. 



Gene Ontology functional enrichment analysis of DEGs in 
the drought-tolerant genotype 

GO (Gene Ontology), an internationally standardized gene 
function classification system, is widely used for gene 
functional annotation and enrichment analysis. GO en- 
compasses three domains: cellular component, biological 
process, and molecular function. The basic GO unit is the 
GO term. Every GO term belongs to a particular category. 
In DGE analysis, GO functional enrichment analysis is 
performed using hypergeometric testing to map all DEGs 
to terms in the GO database (http://www.geneontology. 
org/). To identify significantly enriched GO terms, we cal- 
culated the probability (P) of a particular GO term being 
significantly enriched in a given set of genes according to 
the following formula: 

m=l (M\ ( N-M \ 
p — 1 \ n-i ) 

~ 2^ (N\ 

i=0 \n) 

where N is the total number of genes, n is the number of 
DEGs in N, M is the total number of genes annotated with 
the particular GO term, and m is the number of DEGs in 
M. GO terms with Bonferroni-corrected P-values < 0.05 
were defined as significantly enriched in DEGs. 

In our study, 6,582 DEGs were divided into 51 func- 
tional groups (Figure 6). With respect to cellular compo- 
nents, most DEGs were involved in the categories of 
cell structure (GO: 0005623), intracellular parts (GO: 
0044424), and cytoplasm (GO: 0005737). Under biological 
process, 68.2% of DEGs were found to be involved in 
metabolic processes (GO: 0008152), with 385 DEGs en- 
coding products involved in carbohydrate metabolic pro- 
cesses (GO: 0005975), 318 encoding products involved in 
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Figure 4 Number of DEGs in two genotypes under dehydration and rehydration treatments. (A) Number of up- or down-regulated genes 
in the drought-tolerant genotype. (B) Number of up- or down-regulated genes in the drought-sensitive genotype. 



lipid metabolic processes (GO: 0006629), 253 having 
products involved in amino acid metabolic processes 
(GO: 0006520), and 5 encoding products involved in 
energy reserve metabolic processes (GO: 0006112). 
There were 896 DEGs associated with response to 
stress (GO: 0080134), including response to water 
deprivation (GO: 0009414), salt stress (GO: 0009651), 
temperature stimulus (cold or heat) (GO: 0009266), os- 
motic stress (GO: 0006970), abscisic acid stimulus (GO: 
0009737), and metal ions (GO: 0010038). In addition, 



DEGs were identified in root and leaf systems, includ- 
ing 79 DEGs related to tap and lateral root develop- 
ment (GO: 0048364) and 18 DEGs associated with leaf 
development (GO: 0048366), of which 10 were involved 
in stomatal movement (GO: 0010118). In regard to mo- 
lecular function, many DEGs were found encoding ion 
binding proteins (GO: 0043167) and receptor signaling 
proteins (GO: 0004702), with the majority involved in 
transcription factor (GO: 0003712) activities and pro- 
tein kinase activities (GO: 0004672) (Additional file 10). 
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Figure 5 Venn diagram of DEGs in two genotypes. (A) Venn diagram of DEGs in leaves under dehydration. (B) Venn diagram of DEGs in roots 
under dehydration. (C) Venn diagram of DEGs in leaves under rehydration. (D) Venn diagram of DEGs in roots under rehydration. 
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Figure 6 Histogram of gene ontology classifications. The three GO domains — biological process, cellular components, and molecular 
function — are shown. Left and right y-axes indicate gene category frequencies for each domain and the number of genes in each 
category, respectively. 
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Pathway enrichment analysis of DEGs in the drought- 
tolerant genotype 

Pathway enrichment analysis is an effective method for 
elucidation of DEG biological functions. Pathway-based 
analysis can aid identification of significantly enriched 
metabolic pathways and signal transduction pathways in 
DEGs through comparison with their whole-genome 
backgrounds [29]. The formula used for this calculation 
was essentially identical to that used in the GO analysis, 
with pathways having a Q-value < 0.05 defined as those 
with significantly DEGs. 

KEGG is the major public pathway-related database. 
In our study, DEGs in the drought-tolerant genotype 
were associated with 21 KEGG pathways (Figure 7). 
The pathways with the greatest numbers of unique 
sequences were metabolic pathways (1,712 members) 
(koOHOO), which are large complexes associated with 
metabolic processes such as biosynthesis of secondary 
metabolites and carbohydrate, lipid, and amino acid me- 
tabolism [30-33]. Other pathways, such as plant hor- 
mone signal transduction (573 members) (ko04075), 



plant-pathogen interaction (529 members) (ko04626), 
and protein processing in endoplasmic reticulum (187 
members) (ko04141), were also main enrichment path- 
ways [34-37] (Additional file 11). We believe that these 
pathways may play an important role in plant response 
to dehydration and rehydration conditions. 

GO functional annotation and pathway enrichment 
analysis of DEGs in the drought-tolerant genotype re- 
vealed that the DEGs involved in the most highly 
enriched biological processes and pathways (Additional 
files 10 and 11) primarily encoded transcription factors, 
protein kinases, and other regulatory proteins. These 
DEGs were more likely to determine drought tolerance 
than those in the drought-sensitive genotype. 

DEGs encoding transcription factors in the drought- 
tolerant genotype 

To analyze regulatory mechanisms of drought-responsive 
genes in the drought-tolerant genotype, 37 genes encoding 
transcription factors (TFs) were identified from Jindou21 
leaves and roots (Table 2). These genes could be divided 
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C: Plant hormone signal transduction 

D: Plant-pathogen interaction 

E: Ribosome 

F: Spliceosome 

G: Starch and sucrose metabolism 
H: Phenylpropanoid biosynthesis 
I: RNA transport 

J: Protein processing in endoplasmic 
reticulum 

K: Flavonoid biosynthesis 

L: Ribosome biogenesis in eukaryotes 

M: Purine metabolism 

N: Ubiquitin mediated proteolysis 

O: Amino sugar and nucleotide sugar 

metabolism 

P: mRNA surveillance pathway 

Q: RNA degradation 
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gingerol 

biosynthesis 

U: Pyrimidine metabolism 



Figure 7 Results of KEGG pathway enrichment analyses of DEGs for 21 main KEGG pathways. X- and y-axes represent pathway categories 
and the number of genes in each pathway, respectively. 



into six major groups. The first group contained seven 
auxin response factor (ARF) genes. Two of these genes 
(Glymal4g40540.1 and Glymal3gl7270.2), identified in 
leaves, were suppressed under dehydration. The remaining 
five genes were identified in roots, four of which 
(Glyma07g03840.1, Glyma01g33420.1, Glyma07gl7810.2, 
and Glymal3g29320.2) were induced under dehydration. 
The second group comprised five genes encoding 
ethylene-responsive factors (ERFs), with four of them 
(Glymal7gl5480.1, Glyma05g05 180.1, Glymal3g30990.1, 
and Glymal8g48 740.1) induced under water deficit condi- 
tions. The third group was made up of MYB TF 
family genes. Four (Glymal4g06750.1, GlymaOlgO 1190.1, 
Glymallg05550.1, and Glyma02g00820.1) of the five 
genes were down-regulated under dehydration. The fourth 
group constituted four zinc-finger protein family genes 
(Glyma20g23550.1, Glyma06g04410.1, Glymal3g01290.1, 
and Glymal0g35940.1), two of which were up-regulated 
and two of which were down-regulated. The fifth group 
was composed of ring-finger family genes; one member 
(Glyma20g26780.1) was induced by dehydration in leaves, 
and two (Glyma04g07980.1 and Glymal3g20210.1) were 
suppressed in roots. The sixth group included three HD- 
ZIP family genes (Glyma02g02630.1, Glyma06g20230.1, 
and Glyma09g39360.1) identified only in roots. Under de- 
hydration conditions, two of these were down- regulated 
and one was up-regulated. The remaining TF genes 
encoded heat shock factors (HSFs) (Glymal3g29760.1 and 
Glymallgl4950.1), members of families such as helix-loop- 
helix (bHLH), (Glymal3g39650.1 and Glymal5gl8580.1), 
AP2/EREBP (Glymal4g34590.1 and Glymal0g34760.1), and 
bZIP (Glymal0g01640.1 and Glyma09g38300.1), and other 
types, including NAC domain protein (Glyma20g31210.2) 
and Dof6 (Glymal3g31560.1). 



DEGs encoding protein kinases (PKs) in the tolerant 
genotype 

To better understand signal transduction processes in 
soybean under dehydration, 19 DEGs encoding protein 
kinases were identified in the tolerant genotype (Table 3). 
Of seven DEGs encoding S/T PKs (Glyma20g38980.1, 
Glymal3g43580.2, Glyma05g07580.1, Glymal7gl2620.1, 
Glyma03g29450.1, Glyma09g33120.1, and Glymal5g07820.2), 
three were up-regulated and four were down-regulated under 
dehydration. Three genes (Glyma05g37260.1, Glymal6g03870.1, 
and Glymallg24410.1) encoding CDPKs were down-regulated, 
and one (Glyma03g29450.1) was up-regulated. Simi- 
larly, three genes (Glyma05g02470.1, Glyma05g31 120.1, 
and Glyma20g27580.1) encoding RPKs were suppressed 
and one gene (Glymal7g32780.1) was induced under 
dehydration. Only one gene (Glyma08gl6670.1) encod- 
ing MAPKs was identified in the tolerant genotype. The 
remaining DEGs (Glyma01g00790.1, Glymal6g06940.1, 
and Glymal7g08270.1) encoded three other PKs. 

DEGs encoding other regulatory proteins in the tolerant 
genotype 

Twenty-four genes encoding other drought-responsive 
proteins were also identified in Jindou21 under dehydra- 
tion (Table 4). Among them, one gene (Glyma20g30720.1) 
encoded abscisic stress ripening-like protein and an- 
other (Glymallgl2060.1) encoded salt tolerance-like 
protein. Four genes (Glymal7g37480.1, Glymal4gl 3480.2, 
Glymal4g37440.1, and Glymallg27720.2) encoded two- 
component systems [38]; two of these genes were down- 
regulated in leaves, and two were up-regulated in roots. 
Genes for seven (Glyma20g38610.1, Glyma20g30100.2, 
Glyma20g30100.2, Glymal5gl5010.1, Glyma02g29020.1, 
Glyma20g27660.1, and Glyma07g27370.1) ATP-binding 
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Table 2 Transcription factors responsive to dehydration only in the drought-tolerant genotype 

Transcription factors (TFs) 



Gene Jindou21 leaves (log 2 ratio) Jindou21 roots (log 2 ratio) Annotation 



Auxin-related protein 






Glyma14g40540.1 -4.32 




Auxin response factor(ARF) 


Glyma13g 17270.2 -2.42 




Auxin response factor(ARF) 


Glyma04g43150.1 


-1.56 


Auxin transporter-like protein 3 


Glyma07g03840.1 


1.84 


Auxin-induced protein AUX22 


Glyma01g33420.1 


1.56 


Auxin-induced protein X10A 


Glyma07g17810.2 


1.26 


AUX1 -like auxin influx carrier protein 


Glyma13g29320.2 


2.35 


ARF domain class transcription factor 


Ethylene-responsive proteins 






Glyma17g 15480.1 3.31 




Ethylene transcription factor 1 A 


Glyma05g05180.1 1.83 




ethylene-responsive transcription factor 1 A 


Glyma13g30990.1 1.03 




ethylene-responsive element binding factor 4 


Glyma18g48740.1 


1.81 


ethylene-responsive element binding protein 1 


Glyma15g08360.1 


-2.25 


ethylene responsive transcription factor 12 


MYB transcription factor family 






Glyma14g06750.1 -1.76 




MYB transcription factor MYB75 


Glyma01g01 190.1 -1.45 




GAMYB-binding protein 


Glyma11g05550.1 


-2.05 


MYB transcription factor MYB50 


Glyma15g03920.1 


8.12 


MYB transcription factor MYB82 


Glyma02g00820.1 


-1.89 


GmMYB29A2 


zinc finger protein 






Glyma20g23550.1 1.198624 




zinc finger protein 


Glyma06g04410.1 -1.65 




zinc finger protein 


Glyma13g01290.1 


3.29 


CONSTANS-like zinc finger protein 


Glyma10g35940.1 -7.94 




C2-H2 zinc finger protein 


Ring-H2 protein 






Glyma20g26780.1 1.12 




RING-H 2 finger protein 


Glyma04g07980.1 


-1.03 


ring finger protein 


Glyma13g20210.1 


-1.64 


RING-finger protein 


HD-ZIP family 






Glyma02g02630.1 


1.05 


homeodomain-leucine zipper protein 57 


Glyma06g20230.1 


2.81 


HD-ZIP I protein 


Glyma09g39360.1 


-1.99 


Homeodomain-related 


Heat shock protein 






Glyma13g29760.1 1.50 




heat shock transcription factor 


Glymal 1g14950.1 


1.16 


heat shock cognate protein 1 


AP2/EREBP family 






Glyma14g34590.1 


1.86 


DREB2 


Glymal 0g34760.1 


-4.07 


RAV-like DNA-binding protein 


Helix-loop-helix protein 






Glymal 3g39650.1 


7.82 


helix-loop-helix protein BHLH21 


Glymal 5g1 8580.1 


-1.91 


Helix-loop-helix DNA-binding 
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Table 2 Transcription factors responsive to dehydration only in the drought-tolerant genotype (Continued) 



bZIP transcription factor family 

Glyma10g01640.1 1.70 

Glyma09g38300.1 -7.94 
NAC family 

Glyma20g31210.2 2.15 
Dof6 family 

Glyma13g31560.1 



-1.20 



bZIP transcription factor bZIP59 

leucine zipper-ef-hand containing transmembrane protein 
NAC domain protein 
Dof6 transcription factor 



Note: Genes with FDR < 0.001, |log 2 ratio| >1 and P-value < 0.01 (**) were identified as differentially expressed genes. Significance at 0.01 level is marked with * 



proteins, three cytochrome P450s (Glyma05g09070.1, 
Glymal8g03210.1, and Glyma08g20690.1), and two 
trypsin-inhibitor-like proteins (Glyma09g29310.1 and 
Glyma09g29330.1) were identified and found to be 
down-regulated under dehydration. Increased expression 
was observed for two genes (Glyma06g41 290.1 and 
Glyma06g47780.1) encoding TIR-NBS-LRR-type disease- 
resistance proteins. In addition, four genes encoding other 
types of proteins were identified (Glyma01g04470.2, 
Glyma07g35110.1, Glyma06g35630.1, and Glymal3g07 110.1); 



one was down-regulated in leaves and the others were 
up-regulated. 

Cluster analysis of DEG expression patterns 

Genes with similar expression patterns are usually func- 
tionally correlated. We performed cluster analysis of gene 
expression patterns using Cluster and Java TreeView soft- 
ware [39,40]. The results of hierarchical cluster analysis of 
DEGs under dehydration and rehydration treatments are 
shown in Figure 8. In the drought- tolerant genotype, 48 



Table 3 Protein kinases responsive to dehydration only in the drought-tolerant genotype 

Protein kinase family (PKs) 



Gene 


Jindou21 leaves (log 2 ratio) 


Jindou21 roots (log 2 ratio) 


Annotation 


S/T PKs 








Glyma20g38980.1 


1.12 




serine/threonine protein kinase 


Glyma13g43580.2 


-2.28 




serine/threonine-protein kinase PBS1 


Glyma05g07580.1 


-2.04 




serine/threonine protein kinase 


Glyma17g1 2620.1 




-2.23 


serine/threonine protein kinase 


Glyma03g29450.1 




1.37 


serine-threonine protein kinase 


Glyma09g33120.1 




-1.79 


serine/threonine-protein kinase cx32 


Glyma15g07820.2 




1.41 


serine/threonine protein kinase 


CDPKs 








Glyma05g37260.1 


-1.69 




calcium-dependent protein kinase 


Glyma03g29450.1 


1.37 




calcium-dependent protein kinase isoform 2 


Glyma16g03870.1 




-1.82 


calmodulin-binding receptor-like cytoplasmic kinase 


Glyma11g24410.1 




-1.23 


calmodulin-binding receptor-like cytoplasmic kinase 


RPKs 








Glyma05g02470.1 




-8.23 


receptor-like protein kinase 


Glyma05g31 120.1 




-1.26 


Somatic embryogenesis receptor kinase-like protein 3 


Glyma17g32780.1 


2.84 




stress-induced receptor-like kinase 


Glyma20g27580.1 




-1.51 


receptor-like protein kinase homolog RK20-1 


PKs 








Glyma01g00790.1 




-1.32 


leucine-rich repeat protein kinase 


Glyma16g06940.1 




1.16 


leucine-rich repeat family protein / protein kinase fami 


Glyma17g08270.1 


-1.83 




protein kinase 


MAPKKK 








Glyma08g16670.1 


-1.47 




mitogen-activated kinase kinase kinase alpha 



Genes with FDR < 0.001, |log 2 ratio| > 1 and P-value < 0.01 (**) were identified as differentially expressed genes. Significance at 0.01 level is marked with " : 
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Table 4 Other regulatory proteins responsive to dehydration only in the drought-tolerant genotype 

Other regulatory proteins 



Gene 



Jindou21 leaves (log 2 ratio) 



Jindou21 roots (log 2 ratio) 



Annotation 



Abscisic acid protein 

Glyma20g30720.1 
Salt stress-responsive protein 
Glyma11g1 2060.1 
two-component systems 



-2.02 
-2.74 



2.23 
1.37 
1.02 
10.73 



Glyma17g37480.1 
GlymaHg 13480.2 
Glyma14g37440.1 
Glyma11g27720.2 
ATP binding proteins 
Glyma20g38610.1 
Glyma20g30100.2 
Glyma12g04490.1 
Glyma15g15010.1 
Glyma02g29020.1 
Glyma20g27660.1 
Glyma07g27370.1 
cytochrome P450s 
Glyma05g09070.1 
Glyma18g03210.1 
Glyma08g20690.1 
TIR-NBS-LRRs 
Glyma06g41 290.1 
Glyma06g47780.1 
trypsin-inhibitor-like proteins 
Glyma09g29310.1 
Glyma09g29330.1 
Other type regulated proteins 
Glyma01g04470.2 
Glyma07g35110.1 1.19 
Glyma06g35630.1 -11.44 
Glyma13g07110.1 



-10.41 
-3.44 



7.94 



4.38 



2.19 



7.92 
1.05 



-1.34 
-1.37 
-2.30 



-1.43 

9.07 

-4.08 
-2.39 

1.41 
2.80 



abscisic stress ripening-like protein 
salt tolerance-like protein 
histone H3.2 

aspartate aminotransferase 
asparagine synthetase 2 
asparagine synthetase 2 

ATP-binding cassette transporter 

ATP binding protein 

ATP binding protein 

ATP citrate lyase alpha subunit 

ATP binding protein 

ATP binding protein 

ATP binding protein 

cytochrome P450, putative 

cytochrome P450 monooxygenase CYP90A15 

cytochrome P450 enzyme 

TIR-NBS-LRR type disease resistance protein 
LRR protein-related protein 

pathogen-inducible trypsin-inhibitor-like protein 
pathogen-inducible trypsin-inhibitor-like protein 

pyruvate dehydrogenase alpha subunit 
pyruvate kinase 
tyrosine aminotransferase 
proline dehydrogenase 



Genes with FDR < 0.001, |log 2 ratio| > 1, and P-value < 0.01 (**) were identified as differentially expressed genes. Significance at 0.01 level is marked with ' 



genes were identified between (JW2L vs. JD2L)n(JR0.5 
L vs. JD2L)n(JR2L vs. JD2L)n(JW2R vs. JD2R)n(JR0.5R 
vs. JD2R)n(JR2R vs. JD2R); 35 genes were found be- 
tween (ZW2L vs. ZD2L)n(ZR0.5 L vs. ZD2L)n(ZR2L 
vs. ZD2L)n(ZW2R vs. ZD2R)n(ZR0.5R vs. ZD2R)n 
(ZR2R vs. ZD2R) in the drought-sensitive genotype. 
Under rehydration conditions, some drought-inducible 
genes (Glyma03g29790, Glymallg33760, Glyma08gl2590, 
Glymallgl6000, and Glyma01g43180) were down-regulated, 
while drought-repressed genes (GlymallgQ5550, Glymal5gl5200, 
Glymal3g21350, and Glymal3g35820) were up-regulated upon 



re-watering. These genes may facilitate adaptation to 
drought conditions and promote plant growth recovery 
upon rewatering. 

Candidate DEGs for major roles in response to 
dehydration stress 

The seven most differentially expressed genes (|log 2 ratio | > 
8) (Glymal5g03920, Glyma05g02470, Glymal5gl5010, 
Glyma05g09070, Glyma06g35630, Glyma08gl2590, and 
Glymallgl6000) identified by analysis of DEGs in the 
drought-tolerant genotype and by cluster analysis warrant 
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Figure 8 Hierarchical clustering and heat map presentation of DEGs in two genotypes under dehydration and rehydration. The vertical 
dendrogram indicates the relationship among transcripts across tissues and treatments in the hierarchical cluster analysis. Each column represents 
an experimental condition and each row represents a gene. Expression differences are shown in different colors. Red indicates down-regulation 
and green indicates up-regulation. (A) Jindou21. (B) Zhongdou33. 



further study. Functional annotations of these genes 
are given in Tables 2, 3, and 4. Functional analysis of 
these genes should assist efforts to improve soybean 
drought tolerance. 

Evaluation of housekeeping gene stability and 
verification of DGE results by QRT-PCR 

The software package geNorm was used to evaluate the 
stability of expression of housekeeping genes. UKN2 and 
HDC were ranked the most stable in all samples in our 
experiment, while TUA5 and UBQ10 consistently ranked 
poorly. Optimal numbers of housekeeping genes re- 
quired for RT-PCR data normalization were then deter- 
mined by geNorm (Additional file 12); using this 
program, three housekeeping genes (HDC, UKN2, and 
ACT2/7) were selected to normalize gene expression 
levels (V3/4 = 0.122 < 0.15) (Additional file 13). 

To confirm DGE results, QRT-PCR was conducted on 
eight randomly selected DEGs based on transcriptional 
profile analysis (Figure 9). QRT-PCR results agreed with 
the transcriptional profile data for 96 out of 128 (75%) 



data points. Although specific expression values 
obtained using QRT-PCR were not exactly identical to 
fold changes calculated from expression profiles, both 
methods yielded identical expression trends. QRT-PCR 
results ultimately reflected consistency with the tran- 
scriptional profile data. Sequences of specific primers 
used for QRT-PCR are given in Additional file 3. 

Discussion 

In this study, DGE technology, which allows acquisition 
of more DEGs and related biological information that 
microarray-based techniques [41,42], was used to iden- 
tify DEGs of soybean under dehydration and rehydration 
conditions. This is the first reported attempt to identify 
soybean dehydration-responsive TFs, PKs, and other 
regulatory proteins using both drought-tolerant and 
drought-sensitive genotypes. Although identification of 
DEGs in Arabidopsis, rice, and other plants has been 
reported [43-45], there have been few studies in soybean. 
The acquisition of soybean candidate genes would aid 
the study of soybean drought-tolerant molecular 
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Figure 9 Results of QRT-PCR on eight genes in leaves and roots under dehydration and rehydration treatments. (A) Glyma01g40620.1. 
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mechanisms and, more broadly, drought stress responses 
in dicotyledonous crop plants. 

The drought strategy employed here was similar to a 
previously used air-drying method, where rapid dehydra- 
tion was adopted to simulate drought conditions [46]. In 
our experiments, phenotypes exhibited by soybean in re- 
sponse to drought stress mainly involved leaf wilting and 
rolling; leaf etiolation was not observed. During screen- 
ing for drought- tolerant genotypes, we discovered that 
RWC and REC varied significantly, but chlorophyll con- 
tent and photosynthetic rate did not (data not shown). 
We believe that a decrease in water content due to sto- 
matal closure, the primary symptoms of which are leaf 
wilting and rolling, was caused by the rapid dehydration; 
in contrast, there was no change in chlorophyll content 
and photosynthetic rate in the seedlings in response to 
short-term water deficit. A future search for candidate 
genes related to these physiological indices will incorpor- 
ate long-term drought stress conditions. 

There is much evidence that root and/or leaf traits are 
closely related to mechanisms of drought tolerance [16]. 
Root morphology and development is one of the most 
important factors correlated with mechanisms of 
drought tolerance [47]. The relationship between roots 
and leaves during drought stress is controversial. In one 
study, it was proposed that roots, which are more sensi- 
tive to changes in soil environment, could send drought 



stress signals to leaves [48]. On the other hand, other re- 
searchers have pointed out that leaf growth restraint can 
aid survival during drought by minimizing the evapora- 
tive leaf surface area [49]. In our study, soybean leaves 
and roots were thus collected separately for sequencing. 

Because analysis of DEGs under re-watering is import- 
ant for elucidation of molecular mechanisms during the 
recovery process from dehydration to rehydration [44], 
DEGs were identified during the rehydration process in 
our study in addition to dehydration. We found that 
some drought-inducible genes were repressed under re- 
watering treatment, while the expression of genes 
down-regulated under dehydration was promoted dur- 
ing rehydration. 

Most of the genes differentially expressed during rehy- 
dration could not be annotated. Several genes that were 
functionally annotated encode MYB and homeodomain- 
like TFs and may play central roles in plant response to 
dehydration and growth recovery; however, their func- 
tions require further study. 

GO functional annotation of DEGs 

In addition to their identification, DEGs were annotated 
according to biological process using GO functional 
enrichment analysis. The primary enriched biological 
processes associated with DEGs in the drought-tolerant 
genotype were metabolic processes, response to stress, 
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response to stimulus, transcription, hormone-mediated 
signaling pathways, regulation of gene expression, pro- 
tein modification processes, and plant-type cell wall bio- 
genesis (Additional file 10). A large proportion of 
metabolic processes were related to carbohydrate metab- 
olism, which could provide most of the energy required 
for these pathways under dehydration. Other biological 
process categories enriched in DEGs included aspartate 
family amino acid and Acetyl-CoA metabolic processes 
(Glymal4g37440.1, Glymallg27720.2 and Glymal5gl 5010.1) 
that also contribute to energy release [50]. Proline is an im- 
portant osmolyte, with its synthesis derived from glutamine 
metabolic processes; consequently, genes associated 
with glutamine family amino acid and proline metabolisms 
(Glymal4gl3480.2, Glymallg27720.2, Glymal4g37440.1, 
and Glymal3g07110.1) are necessary for response to osmotic 
stress [51]. Dicarboxylic acid metabolism (Glymal4gl3480.2) 
is also reportedly involved in drought stress [52]. In 
addition, two regulatory proteins (GlymaOlgO 1190.1 
and Glymal5g03920.1) were identified in the tolerant 
genotype. These proteins may protect plants by pro- 
moting the formation of a tolerant morphological 
structure to reduce the impact of adversity [53]. In 
addition to genes associated with metabolic processes 
and morphogenesis, DEGs involved in other processes 
were mainly regulatory proteins, such as transcription 
factors and protein kinases. 

Transcription factors and regulatory proteins are im- 
portant classes of genes that may regulate expression of 
downstream drought-responsive genes [54,55]. Under- 
standing the functions, signaling, and metabolic path- 
ways of TFs, PKs, and other regulatory proteins under 
dehydration will help reveal drought response networks 
in soybean. 

Dehydration stress-inducible TF and PK families in 
soybean 

Transcription factors are known to function in the regu- 
lation of target gene expression. A single transcription 
factor can regulate the expression of multiple genes. 
Transcriptional control of stress-response gene expres- 
sion is a crucial component of plant response to a range 
of environmental stresses [56]. Several transcription 
factors have been reported to be related to drought tol- 
erance in plants. In our study, many DEGs encoding 
TF families were identified in the tolerant genotype 
(Table 2). Among these TFs, MYB TFs, HD-ZIP TFS, 
AP2/EREBP TFs, and zinc finger proteins play a central 
role in drought tolerance [57-60]. Ring finger proteins, 
bHLH proteins, and bZIP TFs regulate the stress- 
responsive ABA signaling pathway [61,62]. HSFs are key 
regulators in the induction of defense systems under 
various stresses [63]. In addition, ARFs and ERFs are 
also involved in stress responses [64,65]. Many members 



of the PK family, which plays a key role in signal trans- 
duction [66,67], were also identified in the tolerant geno- 
type (Table 3). Among these, MAPKs are critical signal 
transduction factors. They can modulate the interaction 
of defense pathways activated by biotic (pathogen infec- 
tion) and abiotic (wounding, drought, salt, and cold) 
factors [68] . MAPKs and S/T PKs can increase multistress 
tolerance [69,70]. CDPKs mainly function in the ABA 
signal pathway and Ca 2+ -mediated regulation [67,71], 
and RPKs enhance abiotic stress in Arabidopsis and 
rice [72,73]. We thus believe that the presence of these 
TFs and PKs detected by DGE profiling indicates that 
various signal molecules act to improve drought toler- 
ance in soybean. 

Other dehydration stress-inducible regulatory proteins in 
soybean 

In addition to TFs and PKs, some regulatory proteins re- 
lated to plant growth, development, and metabolism 
were detected in this study (Table 4). It has been 
reported that ATP synthesis during drought stress is es- 
sential for preserving plastid function [74]. Among our 
DEGs, some ATP-binding proteins identified in the tol- 
erant genotype were down-regulated, which might make 
more ATP available for maintenance of cellular activities 
under water deficit. Cytochrome P450s are involved in 
ABA catabolism in drought stress response [75]. They 
can reduce and eliminate reactive oxygen damage, serve 
as precursors to ABA synthesis, and also participate in 
photosynthesis. In our study, the expression of several 
identified cytochrome P450s was suppressed. Mecha- 
nisms involving two-component gene systems play 
important roles in cellular stress signal transduction 
[22]. A two-component histidine kinase 1 (ATHK1) is 
reported to function as an osmosensor, and phospholip- 
ase C (PLC) is thought to function upstream of the 
drought-responsive DREB2 system [76]. In our study, 
four two-component system genes were identified. Two 
genes separately encoding histone and aspartate amino- 
transferase were down- regulated, and another gene encod- 
ing asparagine synthetase was up-regulated. Regulatory 
factors of two-component systems are induced by 
drought, salt, and osmotic stress [77,78]. These identified 
DEGs thus likely play an important role in improving 
drought tolerance. 

Based on pathway enrichment analysis of the DEGs 
using KEGG, the major public pathway- related database 
[79], genes identified in the tolerant genotype are mainly 
involved in signal transduction and metabolic pathways. 

Signal transduction pathways related to TF and PK 
regulation 

Consistent with our results, a key gene encoding indole- 
3-acetaldehyde oxidase has been identified in tomato 
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under water deficit [45]. In our study, we detected seven 
genes involved in auxin/IAA pathways (Table 2). They 
are downstream genes in the tryptophan metabolic path- 
way, possibly associated with auxin/IAA synthesis in re- 
sponse to water deficit. In addition, two genes encoding 
MYC proteins were found to be involved in the jasmonic 
acid signal pathway. The two genes are located down- 
stream of the a-linolenic acid metabolic pathway. Three 
ethylene-responsive factors, involved in downstream 
regulation of cysteine and methionine metabolism, were 
identified as well. Jasmonic acid and ethylene are factors 
involved in response to senescence stress [80,81]. We 
therefore speculate that these five genes encoding MYC 
proteins and ethylene-responsive factors may regulate 
leaf senescence in response to water deficit. With respect 
to the plant-pathogen interaction pathway, identified 
genes included some encoding MYC proteins and PKs. 
MYC proteins are also involved in the jasmonic acid sig- 
naling pathway, and may be key nodes between the two 
pathways. Protein kinases mainly exist in cell mem- 
branes (FLS 2 and BAKIFKK1) and the cytoplasm 
(CDPK, RPM1, and MEKK), and may possibly regulate 
induction of downstream defense -related genes such as 
WRKYs. The stomatal closure signaling pathway, regu- 
lated by ABA, and the immune pathway, dependent on 
FLS 2 and MEKK1, are interconnected [82,83]. In our 
study, we not only detected proteins related to response 
to ABA, but also identified FLS 2 and MEKK1 proteins 
associated with the plant defense pathway. Consequently, 
we suggest that under dehydration these pathogen- 
interaction proteins may increase plant defense abilities in 
addition to enhancing drought tolerance. 

HSP and regulatory enzyme involvement in 
developmental processes and metabolic pathways 

The results of our DGE analysis indicate that some HSP 
proteins and key enzymes may play roles in regulating soy- 
bean growth and development during dehydration. HSP 
genes are molecular chaperones involved in a variety of 
cellular processes, including protein folding, protein trans- 
port across membranes, regulation of protein degradation, 
and prevention of irreversible protein aggregation. HSP 
over-expression increases drought tolerance in transgenic 
rice and tobacco [84,85]. The drought stress response 
functions of HSP genes identified in our study require fur- 
ther validation. Three protein enzymes identified in our 
study— pyruvate dehydrogenase, pyruvate kinase, and 
tyrosine aminotransferase— may be involved in the gly- 
colysis pathway contributing to the release of energy to 
maintain plant normal metabolism under water deficit. 

Conclusions 

This study has demonstrated the usefulness of the DGE 
approach for identification of DEGs between two soybean 



genotypes under dehydration and rehydration conditions. 
Based on these results, seven candidate dehydration- 
responsive DEGs were selected for further study. Future 
work concentrating on their functions under drought 
stress may lead to a better understanding of genetic mech- 
anisms underlying phenotypic differences between the two 
soybean genotypes under dehydration conditions. Elucida- 
tion of such mechanisms is required for the improvement 
of soybean drought tolerance. 

Additional files 



Additional file 1: List of screened soybean material. Fifty seven 
soybean genotypes from different areas of China were used in 
preliminary screening under different dehydration times, from which two 
drought-tolerant genotypes and two drought-sensitive genotypes were 
selected. 

Additional file 2: Tag preparation principles and steps, including 
cDNA synthesis, enzyme digestion, adaptor ligation, and 
sequencing. 

Additional file 3: Primers used in the QRT-PCR analysis. Primers 
listed were used to amplify eight genes that were randomly selected for 
QRT-PCR assays to confirm DEG reliability and five housekeeping genes 
used to quantify gene expression. 

Additional file 4: Categorization and abundance of tags. The 

number of total tags and distinct tags calculated for different tag 
abundance categories from 28 libraries. 

Additional file 5: Sequencing saturation analysis of 28 libraries. The 

number of detected genes was found to increase as the total number of 
tags increased. 

Additional file 6: Genes responsive to dehydration only in tolerant- 
genotype leaves. DEGs between tolerant and sensitive genotypes were 
obtained by bioinformatics analysis; DEGs expressed under drought only 
in the tolerant-genotype leaves are listed. 

Additional file 7: Genes responsive to dehydration only in tolerant- 
genotype roots. DEGs between tolerant and sensitive genotypes were 
obtained by bioinformatics analysis; DEGs expressed under dehydration 
only in the tolerant genotype roots are listed. 

Additional file 8: Genes responsive to rehydration after 
dehydration only in tolerant-genotype leaves. DEGs between 
tolerant and sensitive genotypes were obtained by bioinformatics 
analysis; DEGs expressed under re-watering only in tolerant-genotype 
leaves are listed. 

Additional file 9: Genes responsive to rehydration after 
dehydration only in tolerant-genotype roots. DEGs between tolerant 
and sensitive genotypes were obtained by bioinformatics analysis; DEGs 
expressed under re-watering only in tolerant-genotype roots are listed. 

Additional file 10: Gene Ontology functional annotation of DEGs in 
the drought- tolerant genotype. The most significantly enriched GO 
biological processes and their associated DEGs are listed. 

Additional file 11: Pathway enrichment analysis of DEGs. Pathway 
enrichment analyses of DEGs are summarized. 

Additional file 12: Evaluation of housekeeping gene expression 
stability using geNorm software. Stabilities of TUA, UBQ, ACT2/7, 
UKN2, and HDC housekeeping genes were analyzed. 

Additional file 13: Standardized calculation of a selected gene. The 

expression level of one randomly chosen gene was compared between 
DGE and QRT-PCR results. 



Abbreviations 

ABA: Abscisic acid; ARFs: Auxin-related factors; CDPKs: Calcium-dependent 
protein kinases; DEGs: Differentially expressed genes; DGE: Digital Gene 
Expression Tag profile; ERFs: Ethylene-related factors; HSFs: Heat shock 



Chen et al. BMC Genomics 2013, 14:687 
http://www.biomedcentral.eom/1 471 -21 64/1 4/687 



Page 1 7 of 1 9 



factors; MAPKs: Mitogen-activated protein kinases; PKs: Protein kinases; QRT- 
PCR: Quantitative real-time PCR; REC: Relative electrical conductivity; 
RPKs: Receptor protein kinases; RWC: Relative water content; S^~ PKs: Serine/ 
threonine protein kinases; TFs: Transcription factors; HSP: Heat shock protein; 
MEKK: MEK kinase; FLS 2 : Flagellin sensing2; RPM1: RESISTANCE TO P. 
SYRINGAE PV MACULICOLA 1. 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

XAZ and WBL conceived the experiments, which LMC designed and 
performed. Data analysis was conducted by LMC and WC; RZ provided the 
soybean seeds used for research. LMC, CW, AHS, ZHS, CJZ, DZQ, ZLY, and 
SLC supplied reagents, materials, and analytical tools. LMC wrote the 
manuscript. All authors read and approved the final manuscript. 

Acknowledgements 

The authors are grateful to the Institute of Crop Sciences of the Chinese 
Academy of Agricultural Sciences (CAAS) and the Economic Crops Research 
Institute of the Shanxi Academy of Agriculture Sciences (SAAS) for supplying 
soybean seeds. This work was funded by the Modern Agro-industry 
Technology Research System and the National Genetically Modified 
Organisms Breeding Major Projects (201 1ZX08004). 

Author details 

^il Crops Research Institute of Chinese Academy of Agriculture Sciences, 
Wuhan 430062, China. 2 Key Laboratory of Oil Crop Biology Ministry of 
Agriculture, Wuhan 430062, China. 3 Key Laboratory of Soybean Biology in the 
Chinese Ministry of Education, Northeast Agricultural University, Harbin 
150030, China. 4 Division of Soybean Breeding and Seed, Soybean Research & 
Development Center, CARS (Key Laboratory of Biology and Genetics & 
Breeding for Soybean in Northeast China, Ministry of Agriculture), Harbin 
150030, China. 

Received: 9 September 2012 Accepted: 25 September 2013 
Published: 6 October 2013 

References 

1 . Boyer JS: Plant productivity and environment. Science 1 982, 21 8:443-448. 

2. Passioura J: The drought environment: physical, biological and 
agricultural perspectives. J Exp Bot 2007, 58:1 13-1 17. 

3. Shinozaki K, Yamaguchi-Shinozaki K: Gene networks involved in drought 
stress response and tolerance. J Exp Bot 2007, 58:221-227. 

4. Cushman JC, Bohnert HJ: Genomic approaches to plant stress tolerance. 
Curr Opin Plant Biol 2000, 3:1 1 7-1 24. 

5. Umezawa T, Fujita M, Fujita Y, Yamaguchi-Shinozaki K, Shinozaki K: 
Engineering drought tolerance in plants: discovering and tailoring genes 
to unlock the future. Curr Opin Plant Biotech 2006, 2:1 13-122. 

6. Shinozaki K, Yamaguchi-Shinozakiy K, Sekiz M: Regulatory network of gene 
expression in the drought and cold stress responses. Curr Opin Plant Biol 
2003, 6:410-417. 

7. Bailly C, Audigier C, Ladonne F, Wagner MH, Coste F, Corbineau F, Come D: 
Changes in oligosaccharide content and antioxidant enzyme activities in 
developing bean seeds as related to acquisition of drying tolerance and 
seed quality. J Exp Bot 2001, 52:701-708. 

8. Breton G, Danyluk J, Charron JBF, Sarhan F: Expression profiling and 
bioinformatic analyses of a novel stress-regulated multispanning 
transmembrane protein family from cereals and Arabidopsis. Plant 
Physiol 2003, 132:64-74. 

9. Wang YC, Jiang J, Zhao X, Liu GF, Yang CP, Zhan LP: A novel LEA gene 
from Tamarix androssowii confers drought tolerance in transgenic 
tobacco. Plant Sci 2006, 171:655-662. 

10. Mardis ER: Next-Generation DNA Sequencing Methods. Annu Rev 
Genomics Hum Genet 2008, 9:387-402. 

1 1 . Asmann YW, Klee EW, Thompson EA, Perez EA, Middha S, Oberg AL, 
Therneau TM, Smith Dl, Poland GA, Wieben ED, Kocher JPA: 3' tag digital 
gene expression profiling of human brain and universal reference RNA 
using lllumina Genome Analyzer. BMC Genomics 2009, 10:531-541. 



12. Feng L, Liu H, Liu Y, Lu ZK, Guo GW, Guo SP, Zheng HW, Gao YN, Cheng SJ, 
Wang J: Power of deep sequencing and agilent microarray for gene 
expression profiling study. Mol Biotechnol 2010, 45:101-1 10. 

13. Sorana Morrissy A, Morin RD, Delaney A, Zeng T, Mcdonald H, Jones S, Zhao 
YJ, Hirst M, Marra MA: Next-generation tag sequencing for cancer gene 
expression profiling. Genome Res 2009, 19:1825-1835. 

14. Hegedus Z, Zakrzewska A, Agoston VC, Ordas A, Racz P, Mink M, Spaink HP, 
Meijer AH: Deep sequencing of the zebrafish transcriptome response to 
mycobacterium infection. Mol Immunol 2009, 46:2918-2930. 

1 5. Qi XH, Xu XW, Lin XJ, Zhang WJ, Chen XH: Identification of differentially 
expressed genes in cucumber (< i > Cucumis sativus</i > L.) root under 
waterlogging stress by digital gene expression profile. Genomics 2012, 
99:160-168. 

16. Manavalan LP, Guttikonda SK, Tran LSP, Nguyen HT: Physiological and 
molecular approaches to improve drought resistance in soybean. Plant 
Cell Physiology 2009, 50:1 260-1 276. 

1 7. Tran LSP, Mochida K: Functional genomics of soybean for improvement of 
productivity in adverse conditions. Funct Integr Genomics 2010, 10:447-462. 

18. Schmutz J, Cannon SB, Schlueter J, Ma J, Mitros T, Nelson W, Hyten DL, 
Song QJ, Thelen J, Cheng JL, Xu D, Hellsten U, May GD, Yu Y, Sakurai T, 
Umezawa T, Bhattacharyya MK, Sandhu D, Valliodan B, Lindquist E, Peto M, 
Grant D, Shu SQ, Goodstein D, Barry K, Futrell-Griggs M, Abernathy B, Du JC, 
Tian ZX, Zhu LC: Genome sequence of the palaeopolyploid soybean. 
A/ature 2010, 463:178-183. 

19. King A, Purcell LC, Brye KR: Differential Wilting among Soybean 
Genotypes in Response to Water Deficit. Crop science 2009, 49:290-298. 

20. Larbi A, Mekliche A: Relative water content (RWC) and leaf senescence as 
screening tools for drought tolerance in wheat. Options 
Mediterraneennes. Serie A, Seminaires Mediterraneens 2004, 60:193-196. 

21 . Liu FH, Guo Y, Gu DM, Xiao G, Chen ZH, Chen SY: Salt tolerance of 
transgenic plants with BADH cDNA. Acta Genetica Sinica 1997, 24:54-58. 

22. Le DT, Nishiyama R, Atanabe Y, Mochida K, Yamaguchi-Shinozaki K, 
Shinozaki K, Tran LSP: Genome-Wide Expression Profiling of Soybean 
Two-Component System Genes in Soybean Root and Shoot Tissues 
under Dehydration Stress. DNA Res 201 1, 1 8(1 ):1 7-29. 

23. Bustin SA: Quantification of mRNA using real-time reverse transcription 
PCR (RT-PCR): trends and problems. J Mol Endocrinol 2002, 29(1):23-39. 

24. Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalisation; 
strategies and considerations. Genes Immun 2005, 6(4):279-284. 

25. Ruan W, Lai M: Actin, a reliable marker of internal control? Clin Chim Acta 
2007, 385(1 -2):1 -5. 

26. Thorrez L, Van Deun K, Tranchevent LC, Van Lommel L, Engelen K, Marchal 
K, Moreau Y, Mechelen IV, Schuit F: Using ribosomal protein genes as 
reference: a tale of caution. Plos One 2008, 3(3):e1854. 

27. Vandesompele J, Preter DK, Pattyn F, Poppr B, Roy NV, Paepe AE, Speleman 
F: Accurate normalization of real-time quantitative RT-PCR date by 
geometric averaging of multiple internal control genes. Genome biology 
2002, 3(7):1-11. 

28. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using 
real-time quantitative PCR and the 2(-Delta Delta C(T) method. Methods 
2001,25:402-408. 

29. Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M: The KEGG resource 
for deciphering the genome. Nucleic acids res 2004, 32:D277. 

30. Sudha G, Ravishankar G: Involvement and interaction of various signaling 
compounds on the plant metabolic events during defense response, 
resistance to stress factors, formation of secondary metabolites and their 
molecular aspects. Plant Cell Tiss Org 2002, 71:181-212. 

31. Mizock BA: Alterations in carbohydrate metabolism during stress: a 
review of the literature. Am J Med 1 995, 98:75-84. 

32. Shi X, Mee-Len C: New roles for acyl-CoA-binding proteins (ACBPs) in 
plant development, stress responses and lipid metabolism. Prog Lipid Res 
2011,50:141-151. 

33. Harding HP, Zhang Y, Zeng H, Novoa I, Lu PD, et al: An integrated stress 
response regulates amino acid metabolism and resistance to oxidative 
stress. Mol cell 2003, 1 1 :61 9-633. 

34. Winkel-Shirley B: Biosynthesis of flavonoids and effects of stress. Curr Opin 
Plant Biol 2002, 5:218-223. 

35. Watanabe S, Nakagawa A, Izumi S, Shimada H, Sakamoto A: A RNA 
interference-mediated suppression of xanthine dehydrogenase reveals 
the role of purine metabolism in drought tolerance in Arabidopsis. 
FEBSIett 2010, 584:1181-1186. 



Chen et al. BMC Genomics 2013, 14:687 
http://www.biomedcentral.eom/1 471 -21 64/1 4/687 



36. Bryan J, Aguilar-Bryan L: Sulfonylurea receptors: ABC transporters that 
regulate ATP-sensitive K + channels. Bba-Biomembranes 1999, 

1461:285-303. 

37. Frugis G, Chua NH: Ubiquitin-mediated proteolysis in plant hormone 
signal transduction. Trends Cell Biol 2002, 12:308-31 1. 

38. Hwang I, Chen HC, Sheen J: Two component signal transduction 
pathways in Arabidopsis. Plant Physiol 2002, 129:500-515. 

39. Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display 
of genome-wide expression patterns. Proc Natl Acad Sci USA] 998, 
95:14863-14868. 

40. Saldanha AJ: Java Treeview-extensible visualization of microarray data. 
BIOINFORMATICS 2004, 20:3246-3248. 

41 . Wang QQ, Liu F, Chen XS, Ma XJ, Zeng HQ, Yang ZM: Transcriptome 
profiling of early developing cotton fiber by deep-sequencing reveals 
significantly differential expression of genes in a fuzzless/lintless mutant. 
Genomics 2010, 6:369-376. 

42. Eveland AL, Satoh-Nagasawa N, Goldshmidt A, Meyer S, Beatty M, Sakai H, 
Ware D, Jackson D: Digital gene expression signatures for maize 
development. Plant Physiol 2010, 154:1024-1039. 

43. Oono Y, Seki M, Nanjo T, Narusaka M, Fujita M, Satoh R, Satou M, Sakurai T, 
Ishida J, Akiyama K, Lida K, Maruyama K, Satoh S, Yamaguchi-Shinozaki K, 
Shinozaki K: Monitoring expression profiles of Arabidopsis gene 
expression during rehydration process after dehydration using ca. 7000 
full-length cDNA microarray. Plant J 2003, 34:868-887. 

44. Degenkolbe T, Do PT, Zuther E, Repsilber D, Walther D, Hincha DK, KOhl Kl: 
Expression profiling of rice cultivars differing in their tolerance to long- 
term drought stress. Plant Mol Biol 2009, 69:133-153. 

45. Gong PJ, Zhang JH, Li HX: Transcription profiles of drought-responsive 
genes in Transcription profiles of drought-responsive genes in 
modulating transcription signal transduction, and biochemical pathways 
in tomato. J Exp Bot 2010, 61:3563-3575. 

46. Le DT, Nishiyama R, Watanabe Y, Mochida K, Yamaguchi-Shinozaki K, 
Shinozaki K, Tran LSP: Genome-wide survey and expression analysis of 
the plant-specific NAC transcription factor family in soybean during 
development and dehydration stress. DNA Res 201 1, 18:263-276. 

47. Sharp RE, Poroyko V, Hejlek LG, Spollen WG, Springer GK, Bohnert HJ, 
Nguyen HT: Root growth maintenance during water deficits: physiology 
to functional genomics. J Exp Bot 2004, 55:2343-2351. 

48. Wang D, Pan YJ, Zhao XQ, Zhu LH, Fu BY, Li ZK: Genome-wide temporal- 
spatial gene expression profiling of drought responsiveness in rice. BMC 
Genomics 2011, 12:149-163. 

49. Achard P, Cheng H, Grauwe LD, Decat J, Schoutteten H, Moritz T, Straeten 
DVD, Peng JR, Harberd NP: Integration of plant responses to 
environmentally activated phytohormonal signals. Science 2006, 31 1:91-94. 

50. Kaiser WM, Kandlbinder A, Stoimenova M, Glaab J: Discrepancy between 
nitrate reduction rates in intact leaves and nitrate reductase activity in 
leaf extracts: what limits nitrate reduction in situ? Planta 2000, 
210:801-807. 

51. Zhang CY, Wang NN, Zhang YH, Feng QZ, Yang CW, Liu B: DNA 
methylation involved in proline accumulation in response to osmotic 
stress in rice (Oryza sativa). Genet Mol Res 2013, 12(2):1 269-1 277. 

52. Wang DL, Ye WW, Wang JJ, Song LY, Fan WL, Cui YP: Construction of SSH 
library and its analyses of cotton drought associated genes under 
drought stress. ACTA AGRONOMICA SINICA 2010, 36(12):2035-2044. 

53. Cho SK, Kim JE, Park JA, Eom TJ, Kim WT: Constitutive expression of abiotic 
stress-inducible hot pepper CaXTH3, which encodes a xyloglucan 
endotrans- glucosylase/hydrolase homolog, improves drought and salt 
tolerance in transgenic Arabidopsis plants. FEBS Lett 2006, 580:3136-3144. 

54. Bray EA: Genes commonly regulated by water-deficit stress in 
Arabidopsis thaliana. J Exp Bot 2004, 55:2331-2341. 

55. Trewavas AJ, Malho R: Signal perception and transduction: the origin of 
the phenotype. Plant Cell 1997, 9:1 181-1 195. 

56. Singh KB, Foley RC, Sanchez LO: Transcription factors in plant defense and 
stress responses. Current Opinion Plant Biology 2002, 5:430-436. 

57. Shin DJ, Moon SJ, Han S, Kim BG, Park SR, Lee SK, Yoon HJ, Lee HE, Kwon 
HB, Baek DW, Yi BY, Byun MO: Expression of StMYB1R-1, a Novel Potato 
Single MYB-Like Domain Transcription Factor, Increases Drought 
Tolerance. Plant Physiol 201 1, 155:421-432. 

58. Sugano S, Kaminaka H, Rybka Z, Catala R, Salinas J, Matsui K, Ohme-Takagi 
M, Takatsuji H: Stress-responsive zinc finger gene ZPT2-3 plays a role in 
drought tolerance in petunia. Plant J 2003, 36:830-841. 



Page 1 8 of 1 9 



59. Yu H, Chen X, Hong YY, Wang Y, Xu P, Ke SD, Liu HY, Zhu JK, Oliver DJ, 
Xiang CB: Activated expression of an Arabidopsis HD-START protein 
confers drought tolerance with improved root system and reduced 
stomatal density. Plant Cell 2008, 20:1 134-1 151. 

60. Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, Yamaguchi-Shinozaki K, 
Shinozaki K: Two transcription factors, DREB1 and DREB2, with an EREBP/ 
AP2 DNA binding domain, separate two cellular signal transduction 
pathways in drought- and low temperature- responsive gene expression, 
respectively, in Arabidopsis. Plant Cell 1998, 10:1391-1406. 

61 . Li HM, Sun JQ, Xu YX, Jiang HL, Wu XY, Li CY: The bHLH type transcription 
factor AtAlB positively regulates ABA response in Arabidopsis. Plant Mol 
Biol 2007, 65:655-665. 

62. Xiang Y, Tang N, Du H, Ye HY, Xiong LZ: Characterization of OsbZIP23 as a 
Key Player of the Basic Leucine Zipper Transcription Factor Family for 
Conferring Abscisic Acid Sensitivity and Salinity and Drought Tolerance 
in Rice. Plant Physiol 2008, 148:1938-1952. 

63. Nishizawa A, Yabuta Y, Yoshida E, Maruta T, Yoshimura K, Shigeoka S: 
Arabidopsis heat shock transcription factor A2 as a key regulator in 
response to several types of environmental stress. Plant J 2006, 
48:535-547. 

64. Seo PJ, Xiang FN, Qiao M, Park JY, Lee YN, Kim SG, Lee YH, Park WJ, Park 
CM: The MYB96 Transcription Factor Mediates Abscisic Acid Signaling 
during Drought Stress Response in Arabidopsis. Plant Physiol 2009, 

151:275-289. 

65. Wang YH, Wan LY, Hang LX, Zhang ZJ, Zhang HW, Quan RD, Zhong SR, 
Huang RF: An ethylene response factor OsWR1 responsive to drought 
stress transcriptionally activates wax synthesis related genes and 
increases wax production in rice. Plant Mol Biol 2012, 78:275-288. 

66. Ludwig AA, Romeis T, Jones JDG: CDPK-mediated signalling pathways: 
specificity and cross-talk. J Exp Bot 2004, 55:181-188. 

67. Zhu SY, Yu XC, Wang XJ, Zhao R, Li Y, Fan RC, Shang Y, Du SY, Wang XF, 
Wang FQ, Xu YH, Zhang XY, Zhang DP: Two Calcium-Dependent Protein 
Kinases, CPK4 and CPK11, Regulate Abscisic Acid Signal Transduction in 
Arabidopsis. Plant Cell 2007, 19:3019-3036. 

68. Xiong LZ, Yang YN: Disease Resistance and Abiotic Stress Tolerance in 
Rice Are Inversely Modulated by an Abscisic Acid-lnducible Mitogen- 
Activated Protein Kinase. Plant Cell 2003, 15:745-759. 

69. Mao XG, Zhang HY, Tian SJ, Chang XP, Jing RL: TaSnRK2.4, an SNFI-type 
serine/threonine protein kinase of wheat (Triticum aestivum L), confers 
enhanced multistress tolerance in Arabidopsis. J Exp Bot 2010, 
61:683-696. 

70. Zhang L, Li YZ, Lu WJ, Meng F, Wu CA, Guo XG: Cotton GhMKK5 affects 
disease resistance, induces HR-like cell death, and reduces the tolerance 
to salt and drought stress in transgenic Nicotiana benthamiana. J Exp Bot 

2012, 10:1-17. 

71 . Zou JJ, Wei FJ, Wang C, Wu JJ, Ratnasekera D, Liu WX, Wu WH: Arabidopsis 
Calcium-Dependent Protein Kinase CPK10 Functions in Abscisic Acid- 
and Ca2 + -Mediated Stomatal Regulation in Response to Drought 
Stress. Plant Physiol 201 0, 1 54:1 232-1 243. 

72. Osakabe Y, Mizuno S, Tanaka H, Maruyama K, Osakabe K, Todaka D, Fujita Y, 
Kobayashi M, Shinozaki K, Yamaguchi-Shinozaki K: Overproduction of the 
Membrane-bound Receptor-like Protein Kinase 1, RPK1, Enhances 
Abiotic Stress Tolerance in Arabidopsis. J Biol Chem 2010, 285:9190-9201. 

73. Ouyang SQ, Liu YF, Liu P, Lei G, He SJ, Ma B, Zhang WK, Zhang JS, Chen SY: 
Receptor-like kinase OsSIKI improves drought and salt stress tolerance 
in rice (Oryza sativa) plants. Plant J 2010, 62:316-329. 

74. Atkin OK, Macherel D: The crucial role of plant mitochondria in 
orchestrating drought tolerance. Ann Bot 2009, 103:581-597. 

75. Narusaka Y, Narusaka M, Seki M, Umezawa T, Ishida J, Nakajima M, Enju A, 
Shinozaki K: Crosstalk in the responses to abiotic and biotic stresses in 
Arabidopsis: analysis of gene expression in cytochrome P450 gene 
superfamily by cDNA microarray. Plant Mol Biol 2004, 55:327-342. 

76. Shinozaki K, Yamaguchi-Shinozaki K: Molecular responses to dehydration 
and low temperature: differences and cross-talk between two stress 
signaling pathways. Curr Opin Plant Biology 2000, 3:21 7-223. 

77. Chen LT, Luo M, Wang YY, Wu KQ: Involvement of Arabidopsis histone 
deacetylase HDA6 in ABA and salt stress response. J Exp Bot 2010, 
12:3345-3353. 

78. Wang HB, Liu DC, Sun JZ, Zhang AM: Asparagine synthetase gene TaASNI 
from wheat is up-regulated by salt stress, osmotic stress and ABA. 

J Plant Physiol 2005, 162:81-89. 



Chen et al. BMC Genomics 2013, 14:687 
http://www.biomedcentral.eom/1 471 -21 64/1 4/687 



Page 1 9 of 1 9 



79. He YH, Fukushige H, Hildebrand D, Gan SS: Evidence Supporting a Role of 
Jasmonic Acid in Arabidopsis Leaf Senescence. Plant physiology 2002, 

128:876-884. 

80. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, 
Kawashima S, Okuda S, Tokimatsu T: KEGG for linking genomes to life and 
the environment. Nucleic Acids Res 2008, 36:D480-D484. 

81. Nie HZ, Zhao CZ, Wu GH, Wu YY, Chen YF, Tang DZ: SR1, a Calmodulin- 
Binding Transcription Factor, Modulates Plant Defense and Ethylene- 
Induced Senescence by Directly Regulating NDR1 and EIN3. Plant 
physiology 2012, 158:1847-1859. 

82. Lefert PS, Robatzek S: Plant Pathogens Trick Guard Cells into Opening the 
Gates. Cell 2006, 126:832-834. 

83. Pitzschke A, Schikora A, Hi rt H: MAPK cascade signalling networks in plant 
defence. Curr Opin Plant Biol 2009, 12:1-6. 

84. Sato Y, Yokoya S: Enhanced tolerance to drought stress in transgenic rice 
plants overexpressing a small heat-shock protein, sHSP17.7. Plant Cell Rep 
2008, 27:329-334. 

85. Alvim FC, Carolino SM, Cascardo JC, Nunes CC, Martinez CA, Otoni WC, 
Fontes EP: Enhanced accumulation of BiP in transgenic plants confers 
tolerance to water stress. Plant Physiol 2001, 126:1042-1054. 



doi:1 0.1 1 86/1 471 -21 64-1 4-687 

Cite this article as: Chen et al.: Genome-wide transcriptional analysis of 
two soybean genotypes under dehydration and rehydration conditions. 

BMC Genomics 2013 14:687. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at (^\ RiftMM i rpntral 

www.biomedcentral.com/submit \^ ™omea centra I 



